Constraining New Physics with D meson decays 
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The standard model (SM) prediction for the leptonic D 3 — > ivg decay is slightly smaller than 
the world average experimental measurement, while the semileptonic SM branching ratios of the 
D° and D+ meson decays are higher than the world average. Using the observed leptonic and 
semileptonic branching ratios for the D meson decays, we performed a global analysis to constrain 
non standard interactions which mediate the cs —¥ Iv transition through generic Wilson coefficients 
that can be later apply to obtain the respective bounds over the relevant parameters for a huge 
variety of models beyond the standard model. In particular, we obtain significant bounds for the 
Two Higgs Doublet Model Type-II and Type III, the Left-Right symmetric model and Minimal 
Supersymmetric Standard Model with explicit R-Parity violation. We found more restrictive or 
similar bounds as previous analysis showing the potential of these type of analysis to constrain new 
physics. Taking into account the reported error on the form factors, we show it is not possible to 
reconcile our constraints for the Two Higgs Doublet Model Type-II at 68% confidence level with 
previous observations. 
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Introduction. In spite of the Standard Model (SM) suc- 
cess, now favored by the probable recent discovery of the 
Higgs boson [l], Q , the search of a more fundamental the- 
ory at an energy scale bigger than the electroweak scale 
is still open. Interestingly, low energy scale experiments 
may shed some light in the search for such fundamental 
theory due to their possibility of getting high statistics 
and hence indirect observables of New Physics (NP). We 
will use D meson decays as an illustration. Contrary 
to B meson physics, charmed hadronic states are in the 
unique mass range of 0(2 GeV), which allows for strong 
non perturbative hadronic physics Q . Moreover, the cal- 
culations for the relevant form factors, which parameter- 
ize all QCD effects within the hadronic state, have been 
improved significantly reaching a surprising precision Q ■ 

At low energies, most of the extensions to the Stan- 
dard Model reduce to an effective four Fermi interaction, 
usually called Non Standard Interaction NSI, that can 
be parameterized by a generic coefficient (Fig. [IJ. For 
the leptonic and semileptonic AC = AS* D meson de- 
cays, the new particle state should couple to the leptons 
and the second generation of quarks, leaving such effec- 
tive interaction. Any kind of intermediate state, such as 
scalars, vectors or even tensors, are allowed. Examples 
arc the Two Higgs Doublet Model Type-II (THDM-II) 
and Type III (THDM-III) @, the Left-Right symmet- 
ric model (LR) and Minimal Supersymmetric Standard 
Model with explicit R-Parity violation (MSSM-X) [1,13 > 
also illustrated in Fig. [1] In this work we make, for 
the first time, a global analysis including all possible non 
standard interactions, over generic coefficients that pa- 
rameterize any NP effects in the cs — > Iv transition. We 
use the semileptonic D° — > K~lv and D + — > K°l + v 
decays and the leptonic decays of the D s meson. In 
the past only the D s leptonic decays were considered Q 




FIG. 1. Generic charged current non standard interaction 
between two quarks and the leptonic sector. Some Feynman 
diagrams for models beyond SM involving the cs — > Iv tran- 
sition involved in D meson decays are shown. 



and the possibility of using semileptonic decays was just 
mentioned in |{|. Furthermore, we use these obtained 
model independent constraints on some particular mod- 
els, namely THDM-II and THDM-III, LR and MSSM-X 
and show that it is enough to consider those coefficients 
to have reliable bounds on specific models. 
Model Independent Analysis. The search of new physics 
effects in the leptonic and semileptonic processes of 
mesons has two sources of uncertainty that can not be 
separated; the non perturbative long-distance forces that 
bind quarks forming hadronic states and the determi- 
nation of the free parameters of the SM, i.e. fermions 
masses and CKM matrix elements. The non-perturbative 
QCD effects are parameterized introducing form factors. 
On the other hand effective couplings that correspond 
to short distance interactions could receive non-standard 
contributions. Hence flavour-changing meson transitions 
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in the SM have at least two scales involved, the elec- 
troweak scale that is responsible of the flavour chang- 
ing and the scale of strong interactions [lfj. When NSI 
are considered, we assume that the new physics energy 
scale is higher than the electroweak scale, thus the oper- 
ator product expansion formalism (OPE) is suitable 
since it allows the separation between long-distance (low 
energy) and short-distance (high energy) interactions. In 
the OPE the degrees of freedom corresponding to higher 
energies scales are integrated out [lH , resulting an effec- 
tive Lagrangian theory where all high energy physics ef- 
fects are parameterized by Wilson's coefficients, namely 
the effective couplings multiplying the operators of the 
Lagrangian. In this spirit, the most general non-standard 
effective Lagrangian for a semilcptonic transition as the 
one illustrated in Fig. [I] is: 

-c NP = J2 c I q ^Jj( q - 1 r I p iq2 )-(D L r I p 2 e), (i) 

I=S,V,T 

where the indexes q\ and q 2 represent down-type and 
up-type quarks respectively, £ is the charged lepton fla- 
vor and v its corresponding neutrino. Pi. 2 represent the 
chiral projectors L = (1 — 7 5 )/2 and R = (1 + 7 5 )/2. 
Here, the current operators L's are determined by the 
Dirac field bilinears, namely: Fg = 1, Ty = 7^ and Tt = 
(«/2)[7 M , 7"]. The dimensionless coefficients C 1 ^^ have 
a clean interpretation: they are a measurement of how 
big can the NSI be as compared to the SM current, 
since they are weighted by the Fermi constant Gf- This 
paramctrization technique enables us to test NSI when 
the experiments reach certain precision, and in partic- 
ular to look for NP effects at low energies. We use this 
technique in a combined analysis of leptonic and semilcp- 
tonic D meson decays to constrain scalar (S) and vecto- 
rial (V) NSI which conserve lepton flavor for a cs — > W 
transition. We shall not take into account the tenso- 
rial contribution in our global analysis, mainly because 
those tensorial NSI are expected to be relatively small 
compared to scalar and vectorial form factors due to the 
lack of tensorial resonance with these quantum numbers 
at D meson scale. Hence, the relevant parameters with 
the above considerations arc 
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For the purpose of this work, we will consider 



C t 

them as real parameters, because the observables we are 
interested in arc more sensitive to the real contribution 
rather than the imaginary part. 

The decay rate of D s — > Ivt including the SM La- 
grangian plus the NSI Lagrangian eq. [TJ is thus given 
by 
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TABLE I. Theoretical and experimental branching ratios 
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Decay 


Theoretical 


Experimental 


1 
1 


U — r I\ e v e 


o.Doy/o. 


fo _|_ n O/l \Q/ 
l^o.OO rt U.U4J /o 


2 


D° -> K~n + v» 


3.559% 


(3.30 ±0.13)% 


3 


D+ -» K°e+u e 


9.325%. 


(8.83 ± 0.22)% 


4 


D+ K°a+^ 


9.099% 


(9.2 ±0.6)% 


5 


D+ -»• T + V T 


5.106% 


(5.43 ±0.31)% 


6 


Dt -> 


5.235 x 10" 3 


(5.90 ± 0.33) x 10" 3 



On the other hand, the semileptonic D° —> K l + vi and 
D + K°£ + V£ decay rates are given by 
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where in the later expression we have defined Gy = 
+ C^ L )/2V2 and G s = {C S J V R + C% v R )/2y/2. 
Other constants involved in Eqs. ([3J [3]) are: Gf the 
Fermi constant, mi,m Cl m s ,rriK,Tnr) 3 ,rnr) the masses of 
the leptons, charm and strange quarks, the Kaon and D 
meson respectively as reported by PDG jl3j]. The trans- 



ferred energy is q 
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2m,F>EK and Ek is 



the final energy of the Kaon meson. Its allowed energy 



m K < E K < (m 2 D 



' K 



m 2 ) /2mF>. In addition to 



those constants, other physical inputs are: 
The CKM element V cs . The direct measurement of this 
element comes precisely from the D meson decay mea- 
surement. Since we are interested in possible signals of 
NP in those decays, we use instead the average value from 
PDG, i.e. V cs = 0.97344. 

Hadronic form factors. For the leptonic decay, the decay 
constant f Ds is defined by {Q\s^f ^c\D s {p)) = ifo.Pii- 
It has been computed with a precision of the order of 
2% by the HPQCD collaboration Q. The reported 
value is f Ds = 248 ± 2.5 MeV. On the other hand, 
less is known about fo(q 2 ), f+{q 2 ), which are defined 



via (K\s^c\D) -- 
with A^ 1 = (m 2 D 



f + (q 2 )( PD + PK - Ar + f Q (q 2 )A^ 
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m K )/(m c — m s )fo(q ). We use the parametrization [1 
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,fo(q 2 ) = 
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(1 - q 2 /P) ' 



where q 2 = q 2 /r 



(4) 

/+(0) has been computed using 
three-flavor lattice QCD calculations. It has been found 
/+(0) = 0.73 ± 0.03 ± 0.07 and a = 0.50 ± 0.04 and 
P = 1.31 ±0.07 0. 

We will ignore all radiative corrections since they are 
expected to be below the 1% [lj|. Once all those ele- 



ments are incorporated, we compute the theoretical de- 
cay rates predicted by the SM T*^^ tv , ^ t D0^ fK + l - v and 
^ t ^+^ / .fc°e+u nxm § Wilson's coefficient to zero. From 
the experimental side, D s leptonic decays have been 
measured by a number of experiments, namely CLEO 
17 and Belle [l8j]. Semileptonic decays have been ob- 
served with an integrated luminosity of SlSpb^ 1 fl9j . 
From those measurements it is possible to extract the 
lifetimes for the mesons. They result to be t d o = 
(410.1 ± 1.5) x 1(T 15 s, t d + = (1040 ± 7) x 10" 15 s, and 
to s = (500 ± 7) x 10 -15 s. It is now possible to compute 
our theoretical branching ratio (BR) just as Bf h = r'"rj, 
i labeling the six different processes involving D decays. 
Our results are listed in table U as well as the experi- 
mental BRs as reported by PDG Note that there 
is a mismatch between the theoretical and experimental 
branching ratios: the leptonic D s — > tvg decay is slightly 
smaller than the world average experimental measure- 
ment, while the semileptonic SM branching ratio of the 
D° and D + meson decays are higher than the world av- 
erage. Let us suppose that the disagreement comes from 
new physics effects, hence we can compute the range of 
the Wilson coefficients to match the theory and the ex- 
periment. In order to do so, we perform a simple \ 2 
analysis, with X 2 = J2 t (^ h ~ B ex P) 2 /SB 2 . Here SB is 
the error, where, in addition to the experimental error, 
we have assumed an extra 4% theoretical error added in 
quadratures to the experimental error for the leptonic de- 
cays, while for the semileptonic decays we have included 
an extra 10% theoretical error added in quadratures as 
well. The main source of errors comes from the form 
factors. The reported error in fo s induces a 2% error 
in the theoretical leptonic branching ratio. That is, we 
are adding the double error in order to be conservative 
in our bounds. On the other hand, the reported error 
in the lattice determination of fo(q 2 ) and f+(q 2 ) leads 
to a 7% error in the theoretical semileptonic branching 
ratio. Again, we have rounded the error to 10% to be 
conservative. The results for the relevant Wilson coeffi- 
cients are listed in Table [Til The first set corresponds to 
universal NSI, that is, flavor independent interactions. In 
this case there are four coefficients: C a £ v , C^^f, C^^f 

S LR 

and C s ^ v . Some models may induce Vector as well as 
scalar NSI at the same time, hence, we report model in- 
dependent constraints by varying all four parameters at 
a time. We also report the bounds obtained by varying 
only one parameter at a time, that is, the most restricted 
bound. Other models may depend on the flavor of the 
lepton involved. Since we have only six B th s, we can 
perform the \ 2 analysis only if we assume scalar NSI or 
vectorial NSI at a time. In each case, for electron NSI, 
we use the channels i = 1,3, for muon i = 2, 4, 6 and for 
the tau, only a fit can be performed with i = 5; channel 
i as shown in Table U Results for both cases, scalar and 
vectorial flavor dependent NSI are listed in Table HI 
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FIG. 2. Allowed regions for tan /3 and the mass of the charged 
Higgs to be consistent with the D meson decays at 68% C.L. 
obtained by using the Wilson coefficient and by a complete \ 2 
analysis of the BRs. The LEP limit on the mass of a charged 
Higgs is also plotted (2f§]. 



Model Dependent Analysis: Let us make explicit the re- 
lation between the Wilson coefficient and some specific 
models. 

Two Higgs doublet model (THDM): It is one of the sim- 
plest and economical extensions of the SM, see [2oT| . [2l| 
for a review. THDM introduces an additional scalar dou- 
blet that induces scalar charged currents (H^), two neu- 
tral scalar fields and a pseudoscalar neutral field (h°,H° 
and j4°). For D meson decays, the only two parameters 
involved are the new scalar mass (m#+) and the ratio 
of the vacuum expectation values tan (3 of the two Higgs 
doublets. At low energies, the Lagrangian for THDM, in 
the Higgs basis for the c harg e scalars and the mass basis 
for fermions, is given by [22( 



- C H ± = V2/vH + [V Uld]Ul (m Ui XP L + m d YP R )d 3 + 
+ m £ Zv L e R } + H.c. (5) 

with X, Y, Z functions of to#+ and tan f3 different for 
different versions of THDM, and the Wilson coefficients 
will be given by 



C 



S.RR 



v: 



mtm c 



ZX . 



c 



S,LR 
cJLv 



2\/2 



v: s ^zy. 



(0) 



In particular, THDM-II has Natural Flavour Conserva- 
tion, namely the suppression of Flavor Changing Neutral 
Currents (FCNC) at tree level, through a Z2 symmetry 

23j . Interesting bounds have been obtained with meson 
decay experiments [lif and recently LEP has reported a 
lower bound on the mass of the charged Higgs of 80 GeV 

25|. For THDM-II, X = cot (3, Y = Z = tan/3. Now we 
will illustrate the effectiveness of our model independent 
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Universal non standard interactions 





4 pars, at-a-time 68% C.L. 


4 pars, at-a-time 90% C.L. 


1 par. at-a-time 68% C.L. 


1 par. at-a-time 90% C.L. 


n V.LL 

s^iS ,RR 
u sc£v 

U ' sclv 


[-0.748, 0.332] 
[-0.876,0.051] 
[-1.17, 1.10] 
[-1.152, 1.130] 


[-0.840, 0.362] 
[-0.969,0.120] 

[-1.24, 1.11] 
[-1.173, 1.178] 


[-0.014, 0.086] 
[-0.142,-0.039] 
[0.003, 0.010] 
[-0.010,2.7 x 10~ 3 ] 


[-0.048,0.118] 
[-0.175, -0.006] 
[1.8 x 10^ 4 , 0.013] 
[-0.013, -2.1 x 10~ 4 ] 



Flavor dependent scalar non standard interactions 





Ctiv' P2 G ^ + 68% C.L. 


C^' P2 e 5R+ 90% C.L. 


C s Jf P2 GK68%C.L. 


CfJ}' P2 G ^ 90% C.L. 


(~iS,RR , (~<S,LR 

fiS,RR 

pS,LR 
(~iS,RR (~tS,LR 

^ ' SCTV T ^ SCTf T 


[0.0, 0.465] 
[3.9 x 10~ 4 , 0.130] 
[0.0,0.147] 
[-8.79 x 10^ 3 , 0.116] 


[0.0, 0.658] 
[0.0,0.134] 
[0.0,0.150] 
[-5.08 x 10^ 2 , 0.154] 


[-0.468, 0.467] 
[-1.23,0.990] 
[-1.23,0.990] 
[-8.79 x 10^ 3 , 0.116] 


[-0.661,0.660] 
[0.312, .234] 
[0.312, .234] 
[-5.08 x 10^ 2 , 0.154] 


Flavor dependent vectorial non standard interactions 




CX 1,P2 e 5R + 68% C.L. 


C^' P2 G &+ 90% C.L. 


cX 1 ^ 2 e ^ 68% C.L. 


cX 1,P2 e 90% C.L. 




[0.0, 0.056] 
[0.0,0.131] 
[0.0, 0.046] 
[-0.014,0.184] 


[0.0,0.111] 
[0.0,0.178] 
[0.0, 0.082] 
[-0.083,0.245] 


[-0.149,0.041] 
[-0.057,0.168] 
[-0.225,3.0 x 10" 3 ] 
[-0.014,0.184] 


[-0.216,0.101] 
[-0.107, .213] 
[-0.270, 0.045] 
[-0.083, 0.245] 



TABLE II. Model independent constraints at 90% and 68% C.L. We have assumed a 10% error of the total semileptonic 
branching ratio and an extra 4% of the leptonic branching ratio, as a theoretical error, which was added in quadratures to the 
experimental error. 



bounds once we apply them to the Wilson coefficients 
of THDM, eqs. |51 There is a flavor dependence com- 
ing from the mass of the leptons involved. Since this is 
an scalar interaction, we can use the bounds on flavor 
dependent scalar NSI. From Cf c '^ - Cf^E we get the 
region -1.8 x 10~ 3 GcV -1 < (m c - m s tan 2 fi)/M% < 
0.023 GcV" 1 at 68% C.L., which gives the outer re- 
gion of an ellipse and the inner region of an hyperbole 
in the plane (m_H-,tan/3) illustrated in Fig. [2l Fur- 
thermore, with the help of the respective bounds for 



3.9 x 10~ 4 < 



2V2V C 



cs Mi 



< 0.130 



C.L.), an upper limit on Mh is obtained. Together with 
the bounds for Cf^ffl the trapezoidal region in Fig. [5] is 
obtained. The combination of these regions are in excel- 
lent agreement with the region obtained by a complete 
X 2 analysis, which is plotted in Fig. [2] in a shadow grey 
area. This agreement illustrates the effectiveness of us- 
ing generic Wilson coefficient to constrain the relevant 
parameters of models beyond the SM. For completeness, 
the THDM-III correspond to the following definitions: 

x = cot(j T^fe 77 ^ ( fl " 22 + ^ f " 21 + ^ 

Y = tan/3 - -^mj 1 (y 2 % + ^Y 2 % + YfcY. 



and Z 



tan/3 2S££_ 



.23 



rd 

1 2,32 



l Y" 



, where Y[ m 

The corre- 



are the 



NSI, we will consider SM's extensions based on ex- 
tending the SM gauge group including a gauge SU (2) r 



In this case, the Wilson coefficients, that 



: -0.048 < C t 
0.118 and -0T75 < 
{Mw/M^f) < -0.006 



can be deduced from 

2 



[321 ] . are 

V*sm^((M w /M^) 2 - 1) < 
-V^W(2 5i )sin(2£)(l- 



V,LL 
sclv 



c 



sctv 



at 90% C.L., which is consistent with TWIST results 

MSSM-j{: R-Parity violating (RPV) superpotentials 
lead to flavor violating interactions in the leptonic and 
hadronic sector that are usually present in many models 
beyond the Standard Model. A vast majority of observ- 
ables have been used to set the corresponding bounds 
to these effective couplings (for a complete review see 
Q and references therein); in particular, for D meson 
decays [35l - l40l ] . We relate the corresponding Wilson co- 
efficients constrained in the global analysis to the RPV 
couplings which constructively interfere with the Stan- 
dard Model, i.e. through the exchange of a —1/3 electri- 
cally charged squark in a t-channcl, which fixes the neu- 
trino flavor. The only non-vanishing Wilson coefficient is 
C^t L = 2V2V cs /G F j: k \K 2k \ 2 /m 2 r Using the conser- 
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12*17™^. < 4-7x10 



VVaGf 2 - 
Yukawa elements as were defined in [2 
sponding bounds obtained via eqs. |5] for THDM-III are 

interesting since they show relations between Y~^ , /? and a g re e with those found in [f 
the mass of the charged Higgs. 

Left-right model: As an example of a pure vectorial 



vative bounds for the model independent constraints (ta- 
ble HH we get the following constraints at 90% confidence 



level: £ fc |A 



7.54X10" 7 , EJA^lVm 2 



-7 



<1V 



< 1.04 x 10 



|A' m l7m^ < 
~ 7 . Our bounds 



Conclusions. We performed a model independent anal- 
ysis of the non standard interactions in the combined pro- 
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cesses D s — s- Ivi and D° — > K£i>£, using the latest exper- 
imental measurements of the branching ratios compared 
to the theoretical expectations. We have found the corre- 
sponding bounds for the Wilson coefficients that param- 
eterize the contribution of new physics as non standard 
interactions. Those constraints can be applied to any 
model of new physics. In particular we have used them 
to constrain the THDM-II, Type III, the Left-Right sym- 
metric model and MSSM-X- We found more restrictive 
or similar bounds as previous analysis, showing the po- 
tential of these type of analysis to constrain new physics. 
Taken into account the reported error on the form fac- 
tors, we found that our results are in disagreement, at 
68% C.L., with previous constraints on the Two Higgs 
Doublet Model Type-II coming from other processes [25| . 
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